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MINIMI Z ATION OF SONIC-BOOM PARAMETERS 
IN REAL AND ISOTHERMAL ATMOSPHERES 

Christine M. Darden 
Langley Research Center 

SUMMARY 

The procedure for sonic -boom minimization introduced by Seebass and George for 
an isothermal atmosphere has been converted for use in the real atmosphere by means 
of the appropriate equations for sonic-boom pressure signature advance, ray -tube area, 
and acoustic impedance. Results of calculations using both atmospheres indicate that 
except for low Mach numbers or high altitudes, the isothermal atmosphere with a scale 
height of 7620 m (25 000 ft) gives a reasonable estimate of the values of overpressure, 
impulse, and characteristic overpressure obtained by using the real atmosphere. The 
results also show that for aircraft design studies, propagation of a known F-function, or 
minimization studies at low supersonic Mach numbers, the isothermal approximation is 
not adequate. 

INTRODUCTION 

The sonic boom is recognized as being one of the major problems confronting the 
advancement of high-speed aircraft. One approach to the reduction of this annoyance is 
through boom -configured aircraft in which low sonic-boom requirements become a funda- 
mental constraint at the outset of aircraft design. This design constraint is embodied in 
the equivalent area distribution which relates the aircraft to its sonic-boom pressure 
signature and therefore specifies the lengthwise development of aircraft lift and volume 
required to produce desired sonic -boom characteristics. Appropriate equivalent area 
distributions can be obtained from the method of Seebass and George (ref. 1), which 
minimizes various parameters of the pressure signature from a cruising aircraft in an 
isothermal atmosphere for a given weight, altitude, and Mach number and defines the 
F-function and equivalent area distribution required to produce this signature. To provide 
the sqme capability in the real atmosphere, the method of reference 1 has been modified 
herein by the appropriate equations (ref. 2) for signature advance, ray-tube area, and 
acoustic impedance. The present study includes pressure signatures in which the impulse, 
initial shock, and maximum overpressure have been minimized. Comparisons of predic- 
tions for these parameters and of pressure signatures, F-functions, and equivalent area 
distributions from both atmospheres are shown. 



SYMBOLS 


Although values are given in both SI and U.S. Customary Units in this report, the 
measurements and calculations for the investigation were made in U.S. Customary Units. 


A 

A e 

Ar, 


B 

C 

D 

F 

G 

H 

h 

I 

K 


ray-tube area 
equivalent area 

equivalent area at equivalent length £ 

speed of sound 

slope of rise in F -function 

height of F -function at balance point 

constant in F-function equation 

Whitham F-function 

area of function to front balance point 

scale height 

airplane altitude 

impulse of pressure signature, \ P 

J P>0 

reflection factor 


dt 


equivalent length of airplane 


M 

P 


Mach number 


pressure perturbation 


2 



Ap 

Ap c 

% 

s 


t 


V 

W 


x,y 


yf 


y r 


z 

a 


overpressure 

characteristic overpressure, 4I/tt 
ambient pressure 
slope of balancing line 
time 

total time between bow and rear shocks in pressure signature 
speed 

airplane weight 
axial distances 

axial position of front area balancing 
position of rear area balancing in F-function 
vertical distance from airplane axis 
advance of acoustic rays 



y ratio of specific heats, 1.4 for air 

r = 2L±i. 

2 

8 Dirac delta function 

X , axial position at which F-function becomes negative 


P 


density 



Subscripts: 


f front shock 

g ground 

h altitude of initial waveform 

J minimum impulse signature 

max maximum pressure in minimum shock signature 

r rear shock 

s minimum shock signature 

so minimum overpressure signature 

A double prime indicates a second derivative with respect to distance. 

METHOD 

At distances from an aircraft where all the characteristics of a pressure signature 
have coalesced into bow and rear shocks of equal strengths, Jones showed that the lower 
bound shock strength required the F -function to have the form of a Dirac delta function 
at x = 0 (ref. 3). He later (ref. 4) extended this conclusion to the bow shocks of N-wave 
signatures in which the bow and rear shocks were of unequal strengths. However, shape 
changes required in current supersonic transports to approach these lower bounds 
resulted in prohibitively high drag penalties (ref. 5). McLean (ref. 6) observed that with 
a sufficiently long airplane, the signature reaching the ground will not have attained its 
N-wave form. Hayes (ref. 7) also pointed out that in the real atmosphere, characteristics 
coalesce more slowly than in the uniform atmosphere and the shape of the signature 
"freezes"; this increases the possibility of the midfield signature shape reaching the 
ground. This phenomenon led George (ref. 8) to consider signatures which differ from the 
basic N-wave form. He showed that the lower bounds for the bow shock of these signa- 
tures also required the F-function to be characterized by a delta function at x = 0. 

The minimizing form of the F-function in these cases and the observation that the 
effective length available for bow shock minimization is reduced when the rear shock is 


4 



also constrained led to the F-function assumed by Seebass and George in their minimiza- 
tion of the entire signature (fig. 1). Most of the equations which follow were developed 
by Seebass and George (ref. 1) and are repeated here for convenience. 

The assumed form of the minimizing F-function is 

F(y) = G6(y) + By + C (O^a) 

F(y) = By - D (X Z y * Q 

In these equations G, B, C, D, and X are unknown coefficients which are deter- 
mined by the cruise conditions of the aircraft, by the prescribed ratio of bow to rear 
shock, and by the parameter of the signature which is being minimized. The Whitham 
function F(y) represents the shape characteristics of the pressure signature and is 
defined in reference 9 in terms of the equivalent area distribution as 




A " 
A e 


(y - x) 


1/2 


dx 


( 2 ) 


Equation (2) is an Abel integral equation which may be inverted to give the function A e (x) 
in terms of the F-function. When this function is evaluated at £, the result is 

A e (£) = 4 J* F(y) (£ - y) 1/2 dy (3) 

Upon substituting the minimizing form of the F-function into equation (3) and integrating, 
the following equation for the development of cross-sectional area is obtained: 


A e (x) = 4GX 1 / 2 + — Bx 5 / 2 + | Cx 3 / 2 - l(x - A)$(C + D)(x - A) 3 / 2 (4) 

e 15 3 3 

where l(x - A) is the Heaviside unit step function. 

If the effects of aircraft wake and engine exhaust are neglected and if the aircraft 
volume is zero at its base, then the area at £ is entirely due to cruise lift or 


•M- _ 4G£ 1 / 2 + M B£ 5 / 2 + - C£ 3 / 2 - - (C + D)(£ - A) 3 / 2 
p y2 15 3 3 


(5) 
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In an isothermal atmosphere = e" h / H , where H is the pressure scale height. If 

P a,g 

this approximation is used in equation (5), then 

_ ;3W _ /3W 

pV 2 yM 2 p yM 2 p e“k/& 

for the isothermal atmosphere. 

Because the linear pressure signal propagates down the ray tube at the local speed 
of sound and each point of the signal advances according to its amplitude, the signal is 
distorted at the ground and would be generally multivalued. Shocks are located in this 
distribution by the area-balancing technique (ref. 10). For the purpose of illustrating 
area balancing, the value of B is zero in figure 2. The first constraint imposed upon 
F(y) is that front area balancing must occur at y = y^, where y^ is the first point at 
which F(y) = C; that is 



For the real atmosphere, the advance (ref. 2) of any point of the signal is given by 

y (2/3) 1 / 2 p a vvy \vy ^ 


where z is the vertical distance of the signal from the airplane axis, z^ is the vertical 
distance of the initial waveform from the aircraft axis, A and are ray-tube areas 
determined from 



( 8 ) 


and r = The initial waveform must be defined away from the airplane because the 

* 2 

F-function can only represent the body shape accurately at several body lengths away and 
the acoustical theory used in describing the propagation of the signal fails near the air- 
craft. The slope of the balancing line is proportional to the reciprocal of the advance at 
any point of the signal; thus 
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s = 


(9) 


(2j3) 1/2 F(y) 



For the rear area balancing to occur between points l and y r then, 

F(y) dy = |(Bj? - D + F(y r j] (y r - £) (10) 


where y r is the unknown second intersection point of the rear area-balancing line with 
F(y). (See fig. 1.) If a cylindrical wake is assumed, F(y) and its integral for y > l 
can be expressed in terms of F(y) for y < l according to reference 11 as follows: 



(y > i) (ii) 


(y > a) (12) 


It is necessary to define F(y) and its integral for y > l in this way for optimization 
problems since aircraft geometry, and thus the Whitham function, can be varied arbi- 
trarily only in the range 05yi{. 

The constraint on the ratio of shocks is given by 


-±- = ^ (13) 

P r D - B£ + F(y r ) 

To insure that y r is an intersection point of F(y) and the balancing line, then 

F (y r ) = s(y r - l) + Bl - D (14) 


and the slope of F(y) at y r must be less than s. 

, Solving the system of equations (5), (6), (10), (13), and (14) provides values for the 
constants G, C, D, A, and y r> Types of signatures studied include flat -topped signa- 
tures in which overpressure is minimized and B = 0 and signatures in which F(y) is 
allowed to rise between y^ and A, with a resulting minimum shock followed by a pres- 
sure rise (fig. 3). The value of B in this form of F(y) may range between 0 and s. 
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With the values of the coefficients known, the minimizing F-function and area distribu- 
tions may be determined by equations (1) and (4). 


To convert F(y) into pressure near the airplane, the following equation is used: 


P\ = yM 2 F (15) 

Pa /h ^z ^ 1 / 2 


The isothermal approximation was used by Seebass and George in the propagation of this 
pressure distribution to the ground. By using linear theory for the real atmosphere, the 
Rayleigh acoustic energy is conserved along acoustic rays and, thus, 



(16) 


Finally, by using a ground reflection factor K, the pressure perturbations are con- 
verted into the ground signature by 

APg = P g K (17) 


DISCUSSION 

Except where indicated otherwise, all results shown are for the following conditions: 
Mach number, 2.7; altitude, 18 288 m (60 000 ft); weight, 272 155 kg (600 000 lb); equiva- 
lent length, 91.44 m (300 ft); reflection factor, 2.0; scale height, 7620 m (25 000 ft); ratio 
of bow -to -rear shock pressure, 1. 

The minimizing F-functions, pressure signatures, and equivalent areas obtained by 
using the same input conditions for the isothermal atmosphere and the real atmosphere 
are shown in figures 4, 5, and 6, respectively. Notice that even though the level of the 
F-function predicted by the isothermal atmosphere is low (fig. 4), the isothermal pressure 
signature (fig. 5) gives an excellent prediction of the shock level. Thus, the error intro- 
duced in using the isothermal atmosphere to calculate the F-function must, in part, be 
canceled out by using the isothermal atmosphere for signature propagation also. The 
lengths of the signatures (fig. 5) do differ; this indicates possible error in those param- 
eters which vary with the length or area of the pressure signature. Use of the isothermal 
substitution for /3 W/pV 2 results in an area distribution whose total area is too low 
to reflect cruise weight requirements at the given Mach number and altitude in the real 
atmosphere. For the same Mach number and altitude then, the isothermal areas shown 
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in figure 6 actually represent a significantly lighter weight airplane. This lower weight 
also explains the lower level of the F-f unction in figure 4. 

The F-functions obtained for both the isothermal and the real atmospheres were 
input into the real atmosphere sonic -boom program of reference 12, and the results from 
this program were compared with the results predicted by the isothermal atmosphere and 
real atmosphere minimization programs. A much lower shock is predicted by the pro- 
gram of reference 12 than by the isothermal minimization program, as shown in figure 7. 
This result supports the statement that the F-function predicted by the isothermal atmos- 
phere represents a smaller airplane, but this F-function, when also propagated through 
an isothermal atmosphere, nearly corrects itself to give good predictions of the shock 
levels for the input conditions. Figure 8 shows that the real atmosphere minimization 
program and the program of reference 12 predict identical pressure signatures on the 
ground. 

The isothermal prediction of shock levels using a scale height of 7620 m (25 000 ft) 
is satisfactory but the resulting equivalent area distributions cannot be used for aircraft 
design studies. The problem arises because the isothermal approximation gives the 
incorrect pressure at flight altitudes, and this leads to an error in total lift. Within the 
isothermal approximation, this error could be alleviated by using the correct pressure 
in the total lift expression (eq. (5)) or by using the scale height which gives the correct 
pressure at that altitude. For the given input conditions, the area distribution can be 
corrected to within 4 percent of the real distribution by using equation (5) and to within 
2 percent by using the scale height, but the resulting shock level predictions are about 
20 percent and 8 percent too large. Thus, it appears that the isothermal atmosphere 
minimization program using a scale height of 7620 m (25 000 ft) is more accurate in pre- 
dicting the shock level. The corresponding area distribution may be improved by the 
factor rea ^A^ isothermal’ ^§ ure 9, this ratio is plotted as a function of altitude 
for a scale height of 7620 m (25 000 ft). After correction, the isothermal area distribution 
differs by less than 8 percent from the real distribution (fig. 10). It should be noted, how- 
ever, that when this corrected area distribution was propagated through the real atmos- 
phere, the resulting shock level was about 5 percent greater than that given by the real 
atmosphere; this emphasizes the sensitivity of shock level to the shape of the equivalent 
area growth. 

For comparison purposes, the far -field minimum impulse results of Jones (ref. 4) 
have also been modified by the appropriate equations for the real atmosphere, and the 
results for the three types of signatures (fig. 3) are shown in figures 11 and 12. Ground 
overpressure as a function of Mach number, equivalent length, altitude, and airplane 
weight is presented in figure 11 for the real and isothermal atmospheres. Values of B 
in the minimum shock signature are taken to be s/2. Smaller values of the shock Ap„ 

fa 
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are obtained in this signature as B — s or as the signature approaches an N-wave. It 
is believed, however, that values of B too close to s give signatures which are indis- 
tinguishable from high-strength N-waves. There is generally good agreement between 
the real atmosphere and isothermal atmosphere results for the shock levels in both the 
minimum overpressure and minimum shock signatures except at the lower Mach numbers 
where the sonic speed gradient in the real atmosphere causes significant changes in the 
ray -tube area (ref. 2) which are not accounted for in the isothermal atmosphere approxi- 
mation. The difference in the predictions for Ap^ , and Ap_ can be explained by 

g,J g,mct* 

the difference in the advance of the two atmospheres. 

The predicted impulses for the real and isothermal atmospheres are presented in 
figure 12. Because of the greater advance, and therefore greater signature length, the 
real atmosphere predictions are higher than the isothermal atmosphere predictions. The 
isothermal predictions for impulse are within 10 percent of the values predicted by the 
real atmosphere with the exception of those at low Mach numbers or high altitudes, where 
significant error occurs. Again, these differences are explained by the greater effect of 
the speed of sound gradient at low Mach numbers and by the error (which increases with 
altitude) in the isothermal advance of the signal. 

In an effort to define a parameter which includes the effects of both shocks and 
impulse, Warren (ref. 13) introduced the characteristic overpressure, defined as 4l/t^. 
where I is the impulse of the pressure signature and is the total time in seconds 
between the bow and rear shocks in the pressure signature. With this definition, the 
characteristic overpressure agrees with the traditional definition of overpressure for the 
far -field N-wave. Figure 13 presents the characteristic overpressures for the minimum 
shock and minimum overpressure signatures. Isothermal atmosphere predictions for 
Ap c give a reasonable estimate of values predicted by the real atmosphere. Real atmos- 
phere predictions for Ap c are generally lower because of the longer signatures in the 
real atmosphere. 


CONCLUDING REMARKS 

The method of Seebass and George for sonic -boom minimization in an isothermal 
atmosphere has been modified for use in a real atmosphere. Comparisons of results for 
both atmospheres indicate that for a given length, weight, altitude, and Mach number, the 
isothermal approximation with a scale height of 7620 m (25 000 ft) gives a good estimate 
of shock strength. The equivalent area distribution resulting from the isothermal atmos- 
phere, however, is unrealistically low but may be improved by a simple factor so that the 
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correct distribution is predicted to within 8 percent. Results also show that for aircraft 
design studies, propagation of a known F-function, or minimization studies at low super- 
sonic Mach numbers, the isothermal approximation is not adequate. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., January 24, 1975. 
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Figure 1.- Assumed form of the F-function in the minimization scheme. 




(a) Minimum overpressure signature. B = 0. 




(c) Minimum impulse signature. 

Figure 3.- Types of signatures considered. Ap c = 4l/tt- 
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(a) F-function corresponding to minimum overpressure signature. 


ft V2 m l/2 



(b) F-function corresponding to minimum shock signature. 


Figure 4.- Minimizing F-functions. M = 2.7; h = 18 288 m (60 000 ft); 
W = 272 155 kg (600 000 lb); j It = 91.44 m (300 ft); P f/ /P r = 1; K = 2 
H = 7620 m (25 000 ft). ' 
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0 100 200 300 ft 

x 

(a) Area distribution corresponding to minimum overpressure signature. 



(b) Area distribution corresponding to minimum shock signature. 

Figure 6.- Minimizing equivalent area distributions. M=2.7; h = 18 288 m 
(60 000 ft); W = 272 155 kg (600 000 lb); SL = 91.44 m (300 ft); 

P f/ P r = K = 2; H = 7620 m (25 000 ft). 
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(b) Minimum shock signature. 


Figure 7.- Predictions obtained with isothermal atmosphere mimimization pro- 
gram and with program of reference 12. M = 2.7; h = 18 288 m (60 000 ft); 
W = 272 155 kg (600 000 lb); t = 91.44 m (300 ft); P f /P r = 1; K = 2; 

H = 7620 m (25 000 ft). 
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(a) Minimum overpressure signature. 
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(b) Minimum shock signature. 

Figure 8.- Predictions obtained with real atmosphere minimization program 
and with program of reference 12. M = 2.7; H = 18 288 m (60 000 ft); 

W = 272 155 kg (600 000 lb); H = 91.44 m (300 ft); P f /P r =1; K = 2; 

H = 7620 m (25 000 ft). 
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Figure 9.- Correction factor for isothermal area distribution when H = 7620 m (25 000 ft). 




(b) Area distribution corresponding to minimum shock signature. 
Figure 10.- Normalized area distributions. 
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(a) Overpressure as a function of Mach number, h = 18 288 m (60 000 ft); 
W = 272 155 kg (600 000 lb); 1 = 91.44 m (300 ft). 


lb/ft 2 N/m 2 
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(b) Overpressure as a function of equivalent airplane length. M = 2.7; 
h = 18 288 m (60 000 ft); W = 272 155 kg (600 000 lb). 

Figure 11.- Comparison of ground overpressures for real and isothermal atmospheres. 

P f /P r =1; K = 2; H = 7620 m (25 000 ft). 
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(a) Impulse as a function of Mach number, h = 18 288 m (60 000 ft); 
W = 272 155 kg (600 000 lb); l = 91.44 m (300 ft). 



(b) Impulse as a function of equivalent length of airplane. M = 2.7; 
h = 18 288 m (60 000 ft); W = 272 155 kg (600 000 lb). 

Figure 12.- Comparison of impulses for real and isothermal atmospheres. 
P f /P r = 1; K = 2; H = 7620 m (25 000 ft). 
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(a) Characteristic overpressure as a function of Mach number, 
h = 18 288 m (60 000 ft); W = 272 155 kg (600 000 lb); 

It = 91.44 m (300 ft). 



(b) Characteristic overpressure as a function of length. M = 2.7; 
h = 18 288 m (60 000 ft); W = 272 155 kg (600 000 lb). 

Figure 13.- Comparison of characteristic overpressures for real and isothermal 
atmospheres. Pf^/P r =1; K = 2; H = 7620 m (25 000 ft). 
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